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We report the preliminary R values for all the 85 energy points scanned in the energy region of 2-5 GeV 
with the upgraded Beijing Spectrometer (BESII) at Beijing Electron Positron Collider (BEPC). On 
average, the uncertainties of the R values we measured are ~ 7%. The new R values has a significant 
impact on the predicted mass of the Higgs (ran) from the global fit to the electroweak data, and will 
also contribute to the interpretation of the E821 g — 2 experiment. 



1 Introduction 

The QED running coupling constant evalu- 
ated at the Z pole, a(M§), and the anoma- 
lous magnetic moment of the muon, = 
(g — 2)/2, are two fundamental imantities to 
test the Standard Model(SM) 0i. a(Mf), 
as one of the three input parameters in the 
global fit to the electroweak data, is sensitive 
to the predicted mass of the Higgs. Theoreti- 
cally, a M is sensitive to large energy scales and 
very high order radiative corrections a. Any 
deviation between the SM predicted value of 
anomalous magnetic moment of the muon, 
o^ M , and that from the experimentally mea- 
sured one, a^ xp , may hint new physics. How- 
ever, the uncertainties in both a(Mg) and 
a^ M are dominated by the hadronic vacuum 
polarization, which cannot be reliably calcu- 
lated but are related to R values through 
dispersion relations i. Here R is the low- 
est order cross section for e + e~ — > 7* — > 
hadrons, which is defined as R = a(e + e~ —> 
hadrons)/cr(e + e~ — > where the de- 

nominator is the lowest-order QED cross sec- 
tion, a(e + e~ — > /i + = a®^ = 47ra 2 /3s. 

Since the uncertainties in a (M§) and 
a^ M are dominated by the errors of the values 
of R in the cm energy range below 5 GeV a, 
it is crucial to significantly reduce the uncer- 
tainties in the R values measured about 20 
years ago with a precision of about 15-20% 
in the energy region of 2-5 GeV au. 



2 R scan with BESII at BEPC 

Following the first R scan with 6 energy 
points in 2.6-5 GeV range done in 1998 □, 
the BES collaboration did a finer R scan 
with 85 energy points in the energy region 
of 2-4.8 GeV. To understand the beam as- 
sociated background, separated beam run- 
ning was done at 24 energy points and sin- 
gle beam running for both e~ and e + was 
done at 7 energy points distributed over the 
whole scanned energy region. Special runs 
were taken at the J / ip to determine the trig- 
ger efficiency. 

The scan was done with BESII, a con- 
ventional collider detector which has been de- 
scribed in detail in ref. 5. 

3 Data Analysis 

The R values from the BESII scan data are 
measured by observing the final hadronic 
events inclusively, i.e. the value of R is deter- 
mined from the number of observed hadronic 
events (N^ s d ) by the relation 

R= N^-N bg -ZiN u -N^ 

where Nb g is the number of beam associated 
background events; Y2i ^Ui Q = e, x) and 
A 77 are the numbers of misidentified lepton- 
pairs from one-photon and two-photon pro- 
cesses events respectively; L is the integrated 
luminosity; 5 is the radiative correction; ehad 
is the detection efficiency for hadronic events. 
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The hadronic event selection identifies 
one photon multi-hadron production from 
all other possible contamination mechanisms. 
Cosmic rays, lepton pair production, two- 
photon process and beam associated pro- 
cesses are the backgrounds involved in our 
measurement. Clear Bhabha events are first 
rejected. Then the hadronic events are se- 
lected based on charged track information. 
Cuts on fiducial, vertex, track fit quality, 
maximum and minimum energy deposition, 
momentum and time-of-flight are applied to 
select hadronic events. Special attention is 
paid to two-prong events. Additional cuts are 
utilized to further reject cosmic raY, Bhabha 
and beam associated backgrounds □. 

The cosmic rays and part of the lepton 
pair production events are directly removed 
by the event selection. The remaining back- 
ground from lepton pair production and two- 
photon processes is then subtracted out sta- 
tistically according to a Monte Carlo simula- 
tion. 

Most of the beam associated background 
events are rejected by a vertex cut. The same 
hadronic event selection criteria are applied 
to the separated-beam data, and the number 
of separated-beam events, N aep , surviving 
these criteria are obtained. The number of 
the beam associated background events, Nb g , 
in the corresponding hadronic event sample 
is given by Nb g = f X N sep , where / is the 
ratio of the product of the pressure at the col- 
lision region times the integrated beam cur- 
rents for colliding beam runs and that for the 
separated beam runs. 

The beam associated background can 
also be subtracted by fitting the event ver- 
tex along the beam direction with a Gaussian 
for real hadronic events and a polynomial of 
degree two for the background. The differ- 
ence between R values obtained using these 
two methods to subtract the beam associated 
background is about (0.3 ~ 2.3)%, depending 
on the energy. This difference was taken into 
account in the systematic uncertainty. 



The integrated luminosity is determined 
by large-angle Bhabha events using only the 
energies deposited in BSC. 

A cross check using only dE / dx informa- 
tion from the MDC to identify electrons was 
generally consistent with the BSC measure- 
ment; the difference was taken into account 
in the overall systematic error of 1.5-2.6%. 

A special effort has been made by the 
Lund group and BES collaboration to de- 
velop the formalism using the basic Lund 
Model Area Law directly for a Monte Carlo 
simulation, which removes the extreme high 
energy approximation in string fragmenta- 
tion in JETSET Q. The final state simulation 
in LUARLW is exclusive as opposed to JET- 
SET's inclusive method, and LUARLW has 
only one free parameter in the fragmentation 
function. Above 3.77 GeV, the production of 
charmed mesons are included in the genera- 
tor based on Eichten Model ill. 

The parameters in LUARLW are tuned 
with R scan data to reproduce distributions 
of kinematic variables such as multiplicity, 
sphericity, angular and momentum distribu- 
tions, etc. 

We find that the same set of parameters 
can be applied to the energy region below 
open charm, and another set of parameters 
can be used for the energies above it. Pa- 
rameters are also tuned point by point for the 
continuum and we find that the detection effi- 
ciencies determined are consistent within 2%, 
which is taken into account in the systematic 
errors. 

Radiative corrections determined using 
four different schemes □ agreed with each 
other to within 1% below charm threshold. 
Above charm threshold, where resonances are 
important, the agreement is within 1-3%. For 
the measurements reported here, we use the 
formalism of Rcf. 9 and include the differ- 
ences with the other schemes in the system- 
atic error of less than 4%. 

The errors from different sources are 
listed in table 1. 
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Table 1. Error Sources for E C m — 3. GeV. Adding the 
systematic and statistic errors in quadrature gives a 
total error of 5.8%. 



Source N had L e had 1 + 5 Stat. 
Err.(%) 3.3 2.3 3.0 1.3 2.5 
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Figure 1. R values below 5 GeV. 



To further improve the measurement of 
R values at BEPC, one needs better perfor- 
mance from the detector and a better handle 
on the uncertainty arising from the hadronic 
event generator, as well as higher machine lu- 
minosity, particularly for the energies below 
3.0 GeV. 

The preliminary R values obtained at 
all 85 energy are graphically displayed in 
Fig. ^, together with the 6 energy points mea- 
sured in the first scan and those measured 
by Markl, 772 and Pluto about twenty years 
ago. The preliminary R values from BESII 
have an average uncertainty of about 7% and 
are slightly lower than that from the previ- 
ously measurements. The two to three factor 
improvement in precision of the R values in 2- 
5 GeV has a significant impact on the global 
fit to the electroweak data for the determina- 
tion of rriH- The preliminary fit results show 
that the predicted itih is significant increased 
with the preferred value lying just above the 
LEP2 excluded region, and the new x 2 pro- 



file of the fit accommodates the_LEP2 bound 
on the mass more comfortably E30. On the 
other hand, BESII R values can also greatly 
contribute to the interpretation of the E821 
g — 2 measurement i. 
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